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Abstract 



We present measurements of the Cabibbo-suppressed decays A+ — > A°K + 
and A+ — > (both first observations), A+ — > S + i^ + 7r _ (seen with large 

statistics for the first time), A+ — > pK + K~ and A+ — > p0 (measured with 
improved accuracy). Improved branching ratio measurements for the decays 
A+ — > E+if+if ~ and A+ — > which are attributed to W-exchange dia- 

grams, are shown. We also present the first evidence for A+ — * H(1690)-fT + 
and set an upper limit on non-resonant A+ — > X + 'K + K~ decay. This analysis 
was performed using 23.6 fb _1 of data collected by the Belle detector at the 
e + e~ asymmetric collider KEKB. 
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I. INTRODUCTION 



Decays of charmed baryons, unlike charmed mesons, are not colour or helicity suppressed, 
allowing us to investigate the contribution of W-exchange diagrams. There are also possible 
interference effects due to the presence of identical quarks. This makes the study of these 
decays a useful tool to test theoretical models that predict exclusive decay rates 0]. 

During the past several years there has been significant progress in the experimental 
study of hadronic decays of charmed baryons. New results on masses, widths, lifetimes and 
asymmetry decay parameters have been published by various experiments 0. However the 
accuracy of branching ratio measurements does not exceed 30% for many Cabibbo-favoured 
modes: for Cabibbo-suppressed and W-exchange dominated decays, the experimental accu- 
racy is even worse. As a result, we are not yet able to conclusively distinguish between the 
decay rate predictions made by different theoretical models. 

In this paper we present a study of A+ baryons produced in the e + e~ — > qq continuum at 
Belle, relying on the excellent particle identification system of the detector to measure decays 
with kaons in the final state. We report the first observation of the Cabibbo-suppressed 
decays A+ — > A°K + and A+ — > Y,°K + , and the first observation of A+ — > H + K + tt~ with 
large statistics. (Here and throughout this paper, the inclusion of charge-conjugate states 
is implied.) We present improved measurements of the Cabibbo-suppressed decays A+ — > 
pK + K~ and A+ — > p<p, and the W-exchange decays A+ — > Yi + K + K~ and A+ — > we 
also report the first evidence for A+ — > H(1690)A' + , and set an upper limit on non-resonant 
A+ — > Y, + K + K~ decay. All branching ratio measurements and limits are preliminary. 

II. DATA AND SELECTION CRITERIA 

The data used for this analysis were taken on the T(4S) resonance and in the nearby 
continuum using the Belle detector at the e + e~ asymmetric collider KEKB. The integrated 
luminosity of the data sample is equal to 23.6 fb _1 . 

Belle is a general purpose detector based on a 1.5 T superconducting solenoid; a detailed 
description can be found elsewhere ||. Tracking is performed by a silicon vertex detector 
(SVD) composed of three concentric layers of double sided silicon strip detectors, and a 
50 layer drift chamber. Particle identification for charged hadrons, important for the mea- 
surement of final states with kaons and/or protons, is based on the combination of energy 
loss measurements (dE/dx) in the drift chamber, time of flight measurements and aerogel 
Cerenkov counter information. For each charged track, measurements from these three sub- 
detectors are combined to form K/tt and p/K likelihood ratios in the range from to 1, 

P(K/tt) = C(K)/(C(K) + £(7r)), P(p/K) = C(p)/(C(p) + C(K)), 

where C(p), C(K) and are the likelihood values assigned to each identification hypoth- 
esis for a given track. 

For the analyses presented here, we require P(K/ir) < 0.9 for pions, P(K/ir) > 0.6 for 
kaons, and P{p/K) > 0.9 for protons unless otherwise stated. Candidate 7r°'s are recon- 
structed from pairs of photons detected in the Csl calorimeter, with a minimum energy of 
50 MeV per photon. Other particles are identified as follows: 
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• A are reconstructed in the decay mode A — > pir , fitting the p and 7r tracks to a 
common vertex and requiring an invariant mass in a ±3MeV/c 2 (~ 3a) interval around 
the nominal value. The likelihood ratio cut on the proton is relaxed to P(p/K) > 0.4. 
We then make the following cuts on the A decay vertex: 

— the difference Z^ st in z-coordinate between the proton and pion at the decay 
vertex must satisfy Zdist < 1 cm; 

— the distance between the vertex position and interaction point (IP) in the r — 
plane must be > 1 mm; 

— the angle a, between the A momentum vector and the vector pointing from the 
IP to the decay vertex, must satisfy cos a > 0.995. 

Finally, to reduce the combinatorial background, we require the A momentum to be 
greater than 1.5 GeV/c. 

• S + are reconstructed in the decay mode S + — > pn°, requiring an invariant mass within 
±10 MeV/c 2 2.0a) of the nominal value. Since the S + mass resolution is dominated 
by the uncertainty in the 7r° momentum, rather than the choice of S + decay vertex, 
we use the IP to estimate the vertex for the 7r° fit; this approximation has been 
checked in the Monte Carlo. The displaced S + vertex is taken into account using the 
impact parameter of the decay proton in the r — plane with respect to the IP, d r $. 
We require the proton to have at least one hit in the SVD, to improve the impact 
parameter resolution, and then make a cut d^p > 500 /zm. 

• S° — > A°7 decays are formed using identified A and photons with calorimeter cluster 
energies > 0.1 GeV; we accept candidates with invariant masses within ±6 MeV/c 2 

1.5a) of the nominal value. To suppress the background from slow photons, the 
cosine of the angle between the 7 in the S° rest frame and the E° boost direction is 
required to be greater than —0.8. 

To suppress combinatorial and BB backgrounds, we require A+ candidates to have scaled 
momentum x p = p* (A^) / \J E^ MS /A — M 2 > 0.5, where p*(A+) and E C ms are the candidate 
momentum and total e + e~ beam energy in the center of mass frame, respectively, and M 
the reconstructed mass of the A+. In modes where there are two or more charged tracks at 
the A+ vertex, we perform a vertex fit and require x 2 / n -d.f. < 9. 

III. FIRST OBSERVATION OF THE DECAYS A+ -» A°K+ AND A+ -» 

The Cabibbo-suppressed decay A+ — > A°K + has not previously been observed. Recon- 
structing A°K + combinations as described in section ||, we see a clear signal at the A+ mass, 
as shown in Fig. |l|(a). 

To study backgrounds due to Cabibbo-allowed decays, we reconstruct a second sample 
with a tight pion identification cut P(K/n) < 0.1 applied to the "kaon", Fig. |l|(b); the 
kaon mass hypothesis is still used in this case. We see a broad structure centered around 
2400 MeV/c 2 : using the MC simulation we find that this is produced by Cabibbo-allowed 
A+ — > A°7r + and A+ — > S°7r + decays. The background due to A+ — > S°7r + is dangerous, 
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M(A°K + ) - M(A°) + 1.1157, GeV/c 2 

FIG. 1. A+ — > A°K + : invariant mass spectra of the selected A°K + combinations. Plot (a) 
is made with the requirement P(K/7r) > 0.6 on the kaon candidate. Plot (b) is made with the 
requirement P(K/tt) < 0.1 to get the form of the contribution from misidentified pions. The other 
selection requirements and the fitting procedure are described in the text. The structure near 
2.4GeV/c 2 is due to misidentified A+ -> A°7r + and A+ -»■ S°7r + decays. 
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since the mass shifts due to the undetected 7 (from E° — > A°7) and 7r — > K misidentification 
partially cancel, leaving a peak close to the A+ mass. Feed-downs from other Cabibbo- 
allowed A c decays, A+ — > A°7r + 7r°, A+ — > A°7r + 7r + 7r~ and A+ — > £°7r + 7r°, produce broad 
reflections due to the missing particle(s), and no prominent structures are seen. We therefore 
fit the distribution using two gaussians with floating central values and widths (to model 
the A+ -> A°7r+ and A+ - 



7T n 



contributions), and a second order polynomial (to model 
the wide reflections and the remaining background); the result of the fit is superimposed on 
Fig. 0(b). All parameters except for the overall normalization are then fixed, and we use 
this function to model the A°7r + background in the main sample (Fig. 0(a)). 

The remaining combinatorial background in Fig. |l](a) is represented using a second 
order polynomial, and the A+ — > A°K + signal is described by a gaussian with width 
a = 5.4MeV/c 2 (fixed from MC); the result of the fit is shown by the superimposed curve. 
We find a yield of 214 ± 30 A+ — > A°K + decays, the first observation of this decay mode. 

For normalization, we use the decay A+ — > A°7r + : the mass distribution of the candidates 
is shown in Fig. |2]. A fit with a gaussian with a floating width for the signal and a second order 
polynomial for the background (restricted to the region above the A+ — ► S°7r + reflection) 
yields 3270 ± 95 events; the fitted signal width o = 6.6 ± 0.2 MeV/c 2 is consistent with 
the MC prediction of 6.8 MeV/c 2 . The relative reconstruction efficiency is found to be 
e(A+ -> A°K + )/e(A+ -»■ A tt + ) = 0.18/0.23 = 0.77 in the Monte Carlo: using this value, we 
extract the branching ratio 



B(A+ - A°K+)/B(Aj 



A 



7T 



0.085 ±0.012 ±0.015. 
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FIG. 2. The invariant mass spectrum for the normalization mode — ► A°7r + . The selection 
requirements and fit are described in the text. The broad feature below the A^" mass is due to 
A+ — > S°7r + decays. 
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The first error is statistical, and the second is systematic, combining the effects of uncertain- 
ties in Kj 7i identification efficiencies and the result of varying the ranges and parameters of 
the fits. 

The Cabibbo-suppressed decay A+ — > TPK + is reconstructed in a similar way, with 
the scaled momentum cut tightened to x v > 0.6 to suppress the large background due to 
soft photons. The invariant mass distribution of the selected TfK + candidates is shown 
in Fig. 3: a peak is seen at the A+ mass, and a reflection due to misidentified two-body 
Cabibbo-allowed A+ decays is seen at higher masses. The superimposed curve shows the 
result of a fit constructed using the method described for A°K + , with the exception that 
the width of the signal gaussian is fixed from the Monte Carlo to a = 5.0 MeV/c 2 in this 
case. We find 70 ± 17 A+ — > TPK + events, the first observation of this decay mode. 

For normalization, we use the decay A+ — > S°7r + , shown in Fig. 4. We fit the dis- 
tribution with a gaussian for the signal, a second gaussian to describe the broad bump 
due to A+ — > A°7r + (with the addition of a random 7), and a second order polyno- 
mial for the remaining background. The mean of the signal gaussian is fixed at the 
A+ mass, and all other parameters are allowed to float: the fitted width of the signal 
a = 5.7 ± 0.2 MeV/c 2 is consistent with the MC prediction a = 6.1 MeV/c 2 . The fit 
gives 1132 ± 39 A+ — > S°7r + decays. The relative reconstruction efficiency is found to be 
e(A+ -> £ fT + )/e(A+ -> S°7r+) = 0.096/0.114 = 0.84 in the Monte Carlo: we then extract 
the branching ratio 

£(A+ -> £ fT + )/i3(A+ -> S°7T + ) = 0.073 ± 0.018 ± 0.016. 

The first error is statistical, and the second is systematic, combining uncertainties from 
particle identification efficiencies and the fitting procedure. 
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M(S°K+) - M(E°) + 1.192, GeV/c 2 

FIG. 3. A+ — > TPK + : invariant mass spectrum of the selected TPK + combinations. The 
selection requirements and fit are described in the text. 
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FIG. 4. The invariant mass spectrum for the normalization mode A+ — ► £°7r + . The selection 
requirements and fit are described in the text. The broad enhancement around the A+ mass region 
is due to A+ — ► A°7r + decays. 



IV. OBSERVATION OF THE A+ -> E+if+Tr" DECAY 



The first evidence for the Cabibbo-suppressed decay A+ — > T l + K + n~ was published by 
the NA32 collaboration in 1992 0: they found 2 events in the signal region. Reconstructing 
Yi + K + ti~ combinations with the cuts of section || tightened to require P(K/n)> 0.9 for the 
kaon, and x p > 0.6, we see a clear signal peak at the A+ mass, as shown in Fig. ^. Tighter 
cuts are used to suppress the large combinatorial background. We also form T l + K + n~ 
combinations using candidates from mass sidebands (> 2er away from the nominal S + 

mass), shown with the shaded histogram: no enhancement is seen near the A+ mass. 



N 

3.0 McV/c 2 




2.22 2.24 2^26 2.28 2^ 2^32 2^4 236 2^38 

Af(S + ir+7T-) - M(£+) + 1.189, GeV/c 2 



FIG. 5. A+ — > T, + K + tt : invariant mass spectrum of the selected Y 1 + K + ir combinations. The 
shaded histogram shows the equivalent spectrum for the S + sidebands. 
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The mass distribution is fitted with a gaussian for the signal (with width fixed to 
3.6 MeV/c 2 from the MC) and a first order polynomial for the background: we find 72 ± 16 
A+ — > T l + K + 7i~ events. For normalization we reconstructed A+ — > £ + 7r + 7r~ decays with the 
same cuts, finding 1432 ± 78 events. The relative efficiency of the A+ — > Y, + K + 7t~ channel 
reconstruction with respect to A+ — > £ + 7r + 7r~ was found to be 0.046/0.055 = 0.85 in the 
Monte Carlo: using this value, we extract the branching ratio 

£(A+ -> Z + K + n-)/B(A+ -> S+7T+7T-) = 0.059 ± 0.014 ± 0.006; 

the systematic error (quoted second) is dominated by the uncertainty in the relative K/ir 
identification efficiency 

V. MEASUREMENT OF THE A+ -» S+A'+A" AND A+ -> DECAYS 

The decays A+ — > Y, + K + K~ and A+ — > proceed dominantly via W-exchange 

diagrams, and were observed by CLEO in 1993 || with branching ratios i3(A+ — ► 
Z + K + K-)/B(At -> j9A"-vr + ) = 0.070 ± 0.012 ± 0.011 and B(A+ -> £+0)/£(A+ -> 
pK~n + ) = 0.069 ± 0.023 ± 0.016. Here we measure these decay channels with improved 
accuracy and provide the first evidence for the A+ — > S(1690)°A' + decay. 

Fig. |^ shows the invariant mass spectrum for A+ — > H + K + K~ combinations selected 
according to section H, with the kaon cuts tightened to P(K/tt) > 0.9, and the impact 
parameter cut for the proton from S + — > pir° relaxed to d r< f, > 200 /mi. A clear peak is 
seen at the A+ mass, over a low background. We fit the distribution using a gaussian (with 
width fixed to 2.2 MeV/c 2 from the MC) plus a second order polynomial: the fit yields 
161 ± 16 A+ — ► H + K + K~ decays. For normalization we reconstructed the A+ — > X + 7r + 7r~ 
decay mode with equivalent cuts, shown in Fig. |7|, and fitted the distribution in the same 
manner: we find 2759 ± 138 A+ — > X + 7r + 7r~ events. The relative efficiency of the A+ — > 
Tj + K + K~ decay reconstruction with respect to the A+ — ■> £ + 7r + 7r~ decay was calculated by 
MC simulation and was found to be 0.052/0.067 = 0.78. We thus extract a branching ratio 

£(A+ -> T l + K + K-)/B(A+ -> S+tt+tt-) = (7.5 ± 0.8 ± 1.5) x 10~ 2 , 

where the second (systematic) error is dominated by the uncertainty in the relative K/ir 
identification efficiency 

In order to obtain the A+ — > X + </> signal, we take H + K + K~ from a ±5 MeV/c 2 win- 
dow around the fitted A+ mass (2286 MeV/c 2 ), and plot the invariant mass of the K + K~ 
combination, Fig. § (points with error bars); the equivalent distribution is also shown for 
Tj + K + K~ in sidebands centred 10 MeV/c 2 below and above the fitted A+ mass (shaded his- 
togram). The distributions are fitted with a Breit-Wigner function (describing the <fi signal) 
convoluted with a gaussian of a fixed width (representing the detector mass resolution) plus 
a 2nd order polynomial multiplied by a square root threshold factor. The intrinsic width of 
the 4> Breit-Wigner function was fixed to its nominal value 0, and the width of the gaussian 
resolution was fixed to 1.0 MeV/c 2 based on the MC simulation. The fit yields 106 ± 12 
events for the </> signal in the A+ region and 15 ± 6 in the A+ sidebands. To extract the 
A+ — > S + contribution we subtract the yield in the sidebands from the yield in the A+ 
signal region, correcting for the phase space factor obtained from the Y, + K + K~ background 



10 



N 

2.0 McV/c 2 



80 
70 
60 
50 
40 
30 
20 
10 




: A 









2.24 



2.26 



2.28 



2.3 



2.32 



2.34 



M{Y, + K + K~) - M(E+) + 1.189, GeV/c 2 
FIG. 6. A+ — > S+K+ivT - : invariant mass spectrum of the selected T l + K + K^ combinations. 
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FIG. 7. The invariant mass spectrum for the normalization mode A+ — ► X + 7r + 7r~, using equiv- 
alent cuts to the A+ — > E+if+i'T - analysis. 
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FIG. 8. Fitting for the A+ — > component: the invariant mass spectra of K + K~ combi- 
nations from the A+ — > Y, + K + K~ signal area (points with error bars) and A+ sidebands (shaded 
histogram) are shown. The selection requirements and fit are described in the text. 



fitting function. After making a further correction for the missing signal outside the A+ 
mass interval, we obtain 93 ± 14 A+ — ► S + decays. 

The relative efficiency of the A+ — > S + </> reconstruction with respect to A+ — ► £ + 7r + 7r~ 
was calculated using the Monte Carlo and found to be 0.050/0.067 = 0.76. Taking into 
account the branching fraction B(<f) — > K + K~) = 49.2% 0, we calculate a branching 
ratio 

£(A+ -> £+</>)/£(A+ -> S+7T+7T-) = (9.1 ± 1.4 ± 1.8) x 10~ 2 . 

The systematic error of 1.8 x 10~ 2 is an estimate based on uncertainties in the relative 
K/tt identification efficiency, and fit variations. In this case there is an additional source of 
systematic error due to the difference in kinematics between the signal and normalization 
modes: this has not yet been taken into account. 

We also searched for the resonant structure in the Yi + K~ system in A+ — > Y, + K + K~ 
decays. Fig. |S] shows the H + K~ invariant mass spectra for Ti + K + K~ combinations 
in a ±5 MeV/c 2 interval around the fitted A+ mass (data points): we also required 
\M(K + K~) — m^l > 10MeV/c 2 to suppress <fi — ► K + K~. Also shown are Y, + K~ from 
Y, + K + K~ combinations selected inside ±5 MeV/c 2 sideband intervals 10 MeV/c 2 below and 
above the fitted A+ mass (shaded histogram). The H + K~ mass distribution shows evi- 
dence for the S(1690)° resonant state. In order to extract this resonant contribution the 
histograms were fitted with a Breit-Wigner function (describing the S(1690)° signal) plus a 
3rd order polynomial multiplied by a square root threshold factor. The fit yields 52.5 ± 15.0 
events for the S(1690)° signal in the A+ region, with a fitted mass and width in good agree- 
ment with previous measurements of the S(1690)° parameters 0. To fit the sidebands, the 
function parameters were fixed to the central values obtained from the signal fit, and the 
normalization was floated: a yield of 7.2 ± 2.8 events was found. 
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FIG. 9. Fitting for the A+ — > H(1690) -fr + component: the invariant mass spectrum of £ + AT~ 
combinations from the A+ — > Ti + K + K~ signal area (points with error bars) and A+ sidebands 
(shaded histogram) are shown, with the (f> — > K + K~ signal region excluded in both cases. The 
selection requirements and fit are described in the text. 

The A+ — > S(1690)°-ft' + contribution was obtained by subtracting the S(1690)° yield in 
the sidebands from the yield in the A+ signal region, correcting the sideband contribution 
using the phase space factor obtained from the Y, + K + K~ background fitting function. After 
a further correction for the missing signal outside the A+ mass interval, we obtained 45 ± 15 
A+ — > S(1690)°-ft" + decays. Assuming the relative efficiency for reconstruction of the A+ — > 
S(1690) i^ + channel with respect to A+ — > £ + 7r + 7r~ to be the same as for the inclusive 
H + K + K~ mode, we find a combined branching ratio 

g( g(A^S+t +l-) +) X ^ 1690 )° ~> = C 2 ' 1 ± °- 7 ± °' 4 ) X 10 ^ 

where we have neglected possible effects due to interference with A+ — > X + </>, etc.. 

Finally, the non-resonant A+ — > Yi + K + K~ contribution is estimated by making invariant 
mass cuts \M(K + K~) - > 10MeV/c 2 and |Af - M H(169 o)| > 20MeV/c 2 to 

suppress the <p and H(1690)° contributions (here, M S ( 1690 ) is the fitted 5(1690)° mass): 
the resulting Y, + K + K~ mass spectrum is shown in Fig. [[(]. A fit with a gaussian (with 
width fixed to 2.2 MeV/c 2 from the MC) plus a first order polynomial, shown on the figure, 
yields 23.8±7.1 events. Integrating the Breit-Wigner function over the allowed M(K + K~) 
region, we found that 14% of the total A+ — > S + signal contributes to this sample: 12.3±2.5 
events. The contribution of the H(1690)° mass tails is estimated to be approximately 12% of 
the fitted S(1690)° signal: 5.4 ± 1.8 events. Subtracting these contributions, 6.1 ± 7.7 non- 
resonant events remain. The phase space correction factor to account for the missing region 
around the <j) and 2(1690)° masses is found to be 1.63 by MC simulation of the non-resonant 
M(K + K~) spectrum: applying this correction we obtain 9.9 ± 12.6 A+ — > H + K + K~ non- 
resonant decays. The relative efficiency for the reconstruction of the A+ — > Y, + K + K~ 
(non- resonant) channel with respect to the A+ — > X + 7r + 7r~ channel was taken to be the 
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FIG. 10. The T, + K + K invariant mass spectrum, after suppressing £ + </> and H(1690)°-KT + 
contributions: the selection requirements and fit are described in the text. 

same as for the inclusive T, + K + K~ mode. Taking into account the systematic error on the 
K/-K identification efficiency, we obtain an upper limit on the branching ratio 

B(At -> ^ + K + K-) non „ res /B(k+ -> S+7T+7T-) < 0.017 

at the 90% confidence level, neglecting possible interference effects between A+ — > S + and 
other resonant decays. 



VI. MEASUREMENT OF THE A+ -> P K + K~ AND A+ -> p<f> DECAYS 

The first evidence for the A+ — > p<j) decay was reported by NA32 in 1990, who claimed a 
signal of 2.8±1.9 events ||. The decay A+ — ► pK + K~ was observed for the first time by E687 
in 1993, who also obtained an upper limit for the branching ratio of A+ — > p<p 0. The most 
recent statistically significant resonant analysis was published by CLEO in 1996, who found 
the following branching ratios: B(A+ -> pK+K~)/B(A+ -> pK~7r + ) = 0.039 ±0.009 ±0.007 
and B(A+ -> p0)/£(A+ -> pK~7T + ) = 0.024 ± 0.006 ± 0.003 g. 

Reconstructing A+ — >■ pK + K~ candidates according to section [T| with the kaon cuts 



tightened to P(K/ir)> 0.9, we see a clear peak at the A+ mass, as shown in Fig. 11. The 
tight P(K/tt) cut is necessary to suppress the large combinatorial background. We fit 
the distribution with a gaussian (with width fixed to 2.8MeV/c 2 from the MC) plus a 
second order polynomial, and find 446 ± 72 A+ — > pK + K~ events. For normalization we 



reconstructed the A+ — > pK~ix + decay mode with equivalent cuts (shown in Fig. 0) and 
fitted the distribution with a double gaussian for the large signal peak, and a second order 
polynomial, finding 33610 ± 1414 events. The relative efficiency of the A+ — > pK~K + decay 
reconstruction with respect to A+ — > pK + 7f~ was found to be 0.21/0.24 = 0.88 in the MC: 
using this value, we extract the branching ratio 

£(A+ -> pK + K-)/B(A+ -> pK~it + ) = (1.50 ± 0.25 ± 0.15) x 10~ 2 . 
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FIG. 11. A+ — > pK + K : invariant mass spectrum of the selected pK + K combinations. 
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FIG. 12. The invariant mass spectrum for the normalization mode — > pif 7r + , using equiv- 
alent cuts to the A+ — > pK + K~ analysis. 
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The systematic error is dominated by the uncertainty in the relative K/tt identification 
efficiency 

In order to obtain the A+ — > p<p signal we take pK + K~ from a ±6MeV/c 2 (±2<t) 
window around the fitted A+ mass (2286 MeV/c 2 ), and plot the invariant mass of the K + K~ 
combination, Fig. |13| (points with error bars); the equivalent distribution is also shown for 
pK + K~ from 6MeV/c 2 sideband intervals 10MeV/c 2 below and above the fitted A+ mass 
(shaded histogram). The distributions are fitted with a Breit-Wigner function (describing 
the signal) convoluted with a gaussian of a fixed width (representing the detector mass 
resolution) plus a 3rd order polynomial multiplied by a square root threshold factor. The 
width of the Breit-Wigner function was fixed to its nominal value 0, and the width of 
the gaussian was fixed to 1.0MeV/c 2 based on MC simulation. The fit yields 430 ± 23 
events for the signal in the A+ region and 232 ± 17 in the sidebands. To extract the 
A+ — > p<p contribution we subtract the yield in the sidebands from the yields in the A+ 
signal region, correcting for the phase space factor obtained from the pK + K~ background 
fitting function. After making a further correction for the signal outside the A+ mass interval 
we obtain 205 ± 30 A+ — > pcj) decays. 
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FIG. 13. Fitting for the A+ — > pep component: the invariant mass spectra of K + K~ com- 
binations from the A+ — > pK + K~ signal area (points with error bars) and sidebands (shaded 
histogram). The selection requirements and fit are described in the text. 

The reconstruction efficiency of the A+ — > pcj) decay relative to A^ — > pK + ir~ was 
calculated using the Monte Carlo and found to be 0.20/0.24 = 0.83. Taking into account 
the branching fraction B{cp — > K + K~) = 49.2% ||, we calculate a branching ratio 

£(A+ -> P (f))/B(At -> pK~ix + ) = (1.50 ± 0.23 ± 0.15) x 10~ 2 . 

The non-0 A+ — > pK + K~ signal is estimated by making an invariant mass cut 
\M(K + K~) — m<f,\ > 10MeV/c 2 to suppress the — > K + K~ contribution: the result- 
ing pK + K~ mass spectrum is shown in Fig [14]. A fit with a gaussian (with width fixed 



to 2.8MeV/c from the MC) and a second order polynomial, shown on the figure, yields 
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FIG. 14. The A+ — > pK + K~ invariant mass spectrum, after suppressing the p0 contribution: 
the selection requirements and fit are described in the text. 

234±59 events. Integrating the Breit-Wigner function over the allowed M(K + K~) region, 
we found that 14% of the total A+ — > p<p signal (27.6 ± 6.2 events) contributes: subtracting 
these, we have 206 ± 59 events. The phase space correction factor accounting for the missing 
region around the mass is found to be 1.075 by MC simulation of the non-0 M(K + K~) 
spectrum: applying this correction we obtain 222 ±64 A+ — > pK + K~ non-0 decays . Assum- 
ing the relative efficiency for reconstruction of A+ — > pK + K~ (non-0) decays with respect 
to A+ — > pK~n + to be the same as for inclusive pK + K~, we calculate a branching ratio 

£(A+ -> pK + K~) non ^/B(Kt -> pK~n + ) = (0.75 ± 0.23 ± 0.08) x 10~ 2 . 

As with previous channels, the uncertainty in the relative Kj n identification efficiency dom- 
inates the systematic error. 

VII. CONCLUSIONS 

In summary, we report the first observation of the Cabibbo-suppressed decays A+ — > 
A°K + and A+ — > H°K + , and the first observation of A+ — > Y, + K + n~ with large statistics. 
The decays A+ — > pK + K~ , A+ — > p<p and A+ — > (pK + K~) m)Q - < j > , and the W-exchange 
decays A+ — > Y, + K + K~ and A+ — > £+0° have been measured with the best accuracy to 
date. We have also observed evidence for the decay A+ — > H(1690)°fT + and set an upper 
limit on the non-resonant decay mode A+ — > Y, + K + K~ . The results for these decay modes 
are listed in Table 1: all values are preliminary. 
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TABLE I. Results obtained in this 



paper. 



Decay mode 



Ratio of branching fractions 



A+ -> A°K+ 
A+ -> 
A+ -» 

A+ -» Y,+K+K- 
A+ -> 

A+ -> ~(1690)°K+ 



(S+K+K-) r 
pK+K- 



A+ 

( P K+^-) non _^ 



c 

A+ 



£(A+ -> A°A^+)/£(A+ -► A°7r+) 

S(A+ -> S°K+)/^(A+ -> S°vr+) 

S(A+ -> S+ J PC+vr-)/i3(A+ -> E+tt+tt") 

B{A+ -> S+ J PC+K-)/,B(A+ -> S+tt+vt-) 

£(A+ -> £+0)/£(A+ -> S+vr+vr-) 

£(A+ -> ~(1690)°K+)x 

x£(~(1690)° -> £+K-)/£(A+ -» S+vr+vr-) 

£(A+ -> Z+K+K-) non _ res /B(A+ -> 

£(A+ -> pK+ J PC-)/i?(A+ -> P K-ir+) 

B{A+ -> p</>)/e(A+ -» pA^+) 

B(A+ -> pK + ir-) non V£(A+ -> pA^+) 



= 0.085 ±0.012 ±0.015 
= 0.073 ±0.018 ±0.016 
= 0.059 ± 0.014 ± 0.006 
= 0.075 ± 0.008 ± 0.015 
= 0.091 ±0.014 ±0.018 

= 0.021 ± 0.007 ± 0.004 
) < 0.017 @ 90% c.l. 
= (1.50 ±0.25 ±0.15) x 10~ 2 
= (1.50 ±0.23 ±0.15) x 10- 2 
= (0.75 ± 0.23 ± 0.08) x 10~ 2 



VIII. ACKNOWLEDGMENT 

We wish to thank the KEKB accelerator group for the excellent operation of the KEKB 
accelerator. We acknowledge support from the Ministry of Education, Culture, Sports, 
Science, and Technology of Japan and the Japan Society for the Promotion of Science; 
the Australian Research Council and the Australian Department of Industry, Science and 
Resources; the Department of Science and Technology of India; the BK21 program of the 
Ministry of Education of Korea and the CHEP SRC program of the Korea Science and 
Engineering Foundation; the Polish State Committee for Scientific Research under contract 
No.2P03B 17017; the Ministry of Science and Technology of Russian Federation; the National 
Science Council and the Ministry of Education of Taiwan; the Japan- Taiwan Cooperative 
Program of the Interchange Association; and the U.S. Department of Energy. 



18 



REFERENCES 



[1] Y. Kohara, Nuovo Cim. Alll, 67-73, 1998, 

M. A. Ivanov et ai, Phys. Rev. D57, 5632-5652, 1998 

K.K. Sharma, R.C. Verma, Phys. Rev. D55, 7067-7074, 1997 

T. Uppal et ai, Phys. Rev. D49, 3417-3425, 1994 
[2] D.E. Groom et al. (Particle Data Group), Eur. Phys. J. C15, 1-878, 2000 
[3] K.Abe et al. (Belle Collaboration), KEK Progress Report 2000-4 (2000), 

to be published in Nucl. Instr. Meth. A 
[4] S. Barlag et al. (NA32 Collaboration), Phys. Lett. B283, 465-470, 1992 
[5] P. Avery et al. (CLEO Collaboration), Phys. Rev. Lett. 71, 2391-2395, 1993 
[6] S. Barlag et al. (NA32 Collaboration), Z. Phys. C 48, 29-45, 1990 
[7] PL. Frabetti et al. (E687 Collaboration), Phys. Lett. B314, 477-481, 1993 
[8] J. Alexander et al. (CLEO Collaboration), Phys. Rev. D53, 1039-1050, 1996 



19 



